'''
Created on 2012-4-23

@author: nzhao
'''
import numpy as np
from math_qm.CoordianteSystem import CoordinateSystem

NV_ZFS                   =  2.87 * 10**6 * 2.0 * np.pi# GHz -> k rad/s
NV_N14_ZFS               = -5.04 * 10**3 * 2.0 * np.pi# MHz -> k rad/s

NV_AXIS = CoordinateSystem(x_axis = np.array([1., -1.,  0.]) * np.sqrt(2), \
                           y_axis = np.array([1.,  1., -2.]) * np.sqrt(6), \
                           z_axis = np.array([1.,  1.,  1.]) * np.sqrt(3)  \
                           )

DIAMOND_LATTICE_CONSTANT = 3.57#Angstrom
DIAMOND_UNIT_CELL_POS    = [ np.array([0,       0,    0]),
                             np.array([1./4, 1./4, 1./4]),
                             np.array([2./4, 2./4,    0]),
                             np.array([3./4, 3./4, 1./4]),
                             np.array([2./4,    0, 2./4]),
                             np.array([0,    2./4, 2./4]),
                             np.array([3./4, 1./4, 3./4]),
                             np.array([1./4, 3./4, 3./4])]